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Abstract 

We prove the existence of a set of two-scale magnetic Wannier orbitals, 
w mn {r), in the infinite plane. The quantum numbers of these states are 
the positions (m, n) of their centers which form a von Neumann lattice. 
Function wqq{t) localized at the origin has a nearly Gaussian shape of 
exp(— r 2 / Al 2 ) / \^2tt for r < v / 2~7r/, where I is the magnetic length. This re- 
gion makes a dominating contribution to the normalization integral. Outside 
this region function, wqo(t) is small, oscillates, and falls off with the Thouless 
critical exponent for magnetic orbitals, r~ 2 . These functions form a conve- 
nient basis for many electron problems. 
71.70.Di, 73.20.Dx 
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I. INTRODUCTION 



Computational methods for interacting electrons in a strong magnetic field have been 
developing rapidly in the last decade.0 They employ different geometries which suggest dif- 
ferent sets of single-electron wave functions and quantum numbers of many-electron states. 
The toroidal geometry was advanced by Yoshioka, Halperin, and Lee.! The symmetry clas- 
sification of the many-electron functions in this geometry was proposed by HaldaneJl After- 
wards most of the finite-size computations were performed in the framework of the spherical 
geometry proposed and developed by Haldane and Rezayi.S 

Landau functions, eigenfunctions of angular momentum, and localized functions called 
coherent stateslH@ are usually used in the plane geometry. Localized functions have obvious 
advantages as applied to the Wigner crystal theory and also to some systems with strongly 
inhomogeneous potential distribution. We expect that these functions may also be effective 
as the basis functions for numerical diagonalization. Since the matrix size increases exponen- 
tially with the number of particles, N, this number is strongly restricted by computational 
facilities, e.g., by iV < 10 for the filling factor v = 1/3. Therefore, it seems tempting to 
develop some procedure for a strong reduction of the dimension of the Hilbert space by 
eliminating functions which make only a small contribution to the low-energy states. We 
believe that the prospects for such an approach are promising since the correlation energy, 
£corr, which can be evaluated by the magnitude of the Laughlin gap or by the energy differ- 
ence between the liquid and solid phases, is small compared with the characteristic Coulomb 
energy, Eq = e 2 /el, where I = (ch/eB) 1 / 2 is the magnetic length. The inequality e COTT <C Eq 
implies that the subspace chosen for diagonalization can be restricted by the electron con- 
figurations with low Hartree-Fock energies. Such an approach has been recently developed 
and successfully applied to electrons on a lattice.! 

As a first step in developing a low-energy states selecting procedure, one should construct 
a basis of localized, particle-like functions with a characteristic size about I. However, some 
general requirements exclude the existence of a complete system of single-electron states 
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which are both exponentially localized and linearly independent. This point should be 
explained in more detail. 

The coherent states^ on a von Neumann lattice are the most localized eigenf unctions for 
an electron in a magnetic field. For such a lattice, with a single flux quantum per a unit 
cell, the area of a unit cell equals 2tc1 2 . The set of coherent states, c mn (r), of the lowest 
Landau level (LLL) can be obtained by the magnetic translations of the function coo(V) = 
exp(— r 2 /4Z 2 )/\/2tF from the origin to all lattice sites, (m, n). The set c mn (r) is complete, 
i.e., an arbitrary function belonging to the LLL can be expanded in functions c mn (r). This 
statement is physically appealing and is supported by a rigorous mathematical proof. 
The functions c mn (r) are obviously non-orthogonal. A straightforward way to orthogonalize 
such a set is to transform it into the Bloch representation, and then go back to the site 
representation. Wannier applied this procedure to a one-dimensional chain of Gaussians.il! 
If the functions c mn (r) were linearly independent, this orthogonalization procedure would 
result in a set of exponentially localized functions. However, Perelomov0 has established 
a non-trivial fact that the set c mn (r) is actually overcomplete by exactly one function and 
presented the explicit form of the linear equation relating these functions; see Eq. (|J) below. 
The overcompleteness of the set c mn (r) imposes hard restrictions on the localization of 
Wannier functions. Zak and collaborators^'!^ have shown that the exponential localization 
and the orthogonality of magnetic orbitals are incompatible. Thoulessli! proved even a 
stronger result by relating the localization of magnetic orbitals to the existence of the Hall 
current. He showed that in systems supporting a Hall current the orbitals should fall off 
with distance no faster than by the inverse-square law, hence, r~ 2 is a critical exponent. 

Having in mind these restrictions on the falloff of magnetic orbitals, one can question 
whether particle-like orbitals suitable for restricted-basis finite-size calculations do exist. It 
is the main result of this paper that a set of two-scale orbitals w mn (r), which satisfy the 
above requirements exists.0 In the small r region, r < \^2ttI, the orbital Woo(r) is very 
close to the Gaussian coo(V), an d this region contributes about 95% to the normalization 
integral. For r > \f2irl, the function u>oo(r) f &ns °ff as r_2 i n exact correspondence with 



the Thouless criterion and shows an oscillatory behavior which ensures orthogonality to the 
different orbitals w mn (r). 

The term Wannier functions is usually applied only to the orbitals which fall off expo- 
nentially Since the functions w mn (r) are normalizable and have a well localized Gaussian 
core, we shall term them as Wannier functions in what follows. 

The paper is organized as follows. In Sec. |T| general equations are presented and the 
behavior of magnetic Bloch functions near the singular point in the momentum plane is 
studied. It is of importance for establishing the completeness of the sets of Bloch and 
Wannier functions. The basic results are discussed in Sec. [II] where two-scale Wannier 
functions w mn (r) of the infinite plane are derived and studied analytically and numerically. 



In Sec. [V| some results for Wannier functions on the finite plaquettes are presented. The 
results are summarized in Conclusion. 



II. GENERAL EQUATIONS 

Let us choose a rectangular unit cell with sides a and b such that ab = 2tt and a 
normalization plaquette with sides L x = aa and L y = (5b, where a and (3 are integers. Here 
and below 1 = 1. It is convenient to define the function c mn (r) centered at a site (m, n) as 

Cmn(r) = T ma T nh c 00 (r) , c 00 (r) = exp(-r 2 /4)/v / 27T , (1) 

where T ma and T nb are operators of magnetic translations!^ along a and b axes, and c mn (r) 
are the eigenfunctions of the Schroedinger equation of an electron in a magnetic field. The 
axial gauge A(r) = ~ z x r is used. The general definition of the operator Tr for a particle 
with a charge e > is: 

T R ^(r) = exp{-^ R ■ A(R)} exp{*r ■ A(R)} V(r - R) • (2) 

Therefore, the explicit form of the function c mn (r) is: 

Cmn(r) = {-) mn exp{-^(r - R mn ) 2 + ~ z • (R mn x r)}/V^ , (3) 
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where R mn = ma + nb. 

The set of functions c mn (r) obeys the Perelomov overcompleteness equation^ 

oo 

£ n m+n c mn (r) = (4) 



mn=-oo 



Therefore functions c mn (r) are linearly dependent. 

The operators T ma T n b in Eq. (P are considered as elements of the outer Kronecker 
product G a (g) Gb of the groups G a and Gb of the magnetic translations along a and b. 
It follows from Eq. (|) that T a T h = T h T a e mb = T h T a . Therefore, T ma and T nh commute 
for arbitrary values of m and n. By definitionjl! the outer product consists of operators 
T m T„b and does not include operators T ma+n b with mn ^ which do not commute with 
the operators of the T ma and T nb types. Therefore, the Kronecker product G a ® Gb forms 
an Abelian group of translations, and one can introduce a two-dimensional momentum k 
and follow the Wannier procedure]^ 

Transformation of the set c mn (r) to the momentum representation results in Bloch func- 
tions: 

oo 

^k(r) = ^2 c mn(r) exp (tkK mn )/^/a(3 z/(k) . (5) 

mn=-oo 

Functions \I/k(r) with different values of k are orthogonal since they belong to different 
irreducible representations of the group G a ®Gb- Values of k can be found from the boundary 
conditions for \l/k( r )- The latter can be chosen either as regular periodic boundary conditions 
with respect to the magnetic translations with periods L x and L y , or from 0-periodic 
conditions when some twist (p = (<f) x , 4> y ) is added.0 Therefore, the components of k can 
take values 



k x = (2ns x + <f) x )/L x , k y = (27cs y + <py)/L y (6) 

inside the Brillouin zone. Here components s x and s y of the vector s = (s x , s y ) are integers, 
s x , s y = 0, ± 1, ... , taking a and (3 values, respectively. Normalization coefficient z/(k) is 
determined by the equation 



z/(k) = a/2tt ^ c mn (0,0) cos(kR m „) . (7) 

mn=-oo 

The inverse Fourier transformation 

W mn (r) = J^*k(r) exp(-ikR mn )/v^ (8) 

k 

results in a set of orthonormal Wannier-type functions: 

oo 

W mn (r) = ^ K^m' - m, n - n) c m / n /(r) , (9) 

m'n'=— oo 

where 

K(f,(m,n) = — y^exp(ikR mra )/A/^(k) . (10) 
a/} k 

The functions K^m, n) and W TOn (r) obey the equations K^ijn + Ma,n + N/5) = 
K^(m, n) exp{z0(M, N)} , and 

(r) = W ron (r) exp{-^(M, N)} , (11) 

where 0(M, N) = M^^. + N0 y , and M and N are integers. Therefore, |W mn (r)| is periodic 
with periods L x and L y for any value of the flux </>. 

The above equation for the Bloch functions ^(r) is known both in LandauiHl and co- 
herent functional basises, and the equation for W mn (r) was discussed more recently as applied 
to finite plaquettes.i These equations have their analogs in the mixed kq representation.0 

Functions \I/k( r ) are orthogonal and, therefore, linearly independent. If they exist for 
any momenta k = (k x , k y ), defined by Eq. (|6]), then the set ^k(r) is complete. Hence, the set 
W mn (r) is complete, too. However, the properties of one of the functions ^k( r ), namely the 
function \Pk (r), where k is the corner of the Brillouin zone, need a special deliberation. At 
this point k x a = ± k y b — ± 7T, and the exponent in Eq. (|5]) equals (— ) m+ ™. Therefore, the 
numerator of \I/ko( r ) turns into zero because of Eq. (f|). If one puts ^k ( r ) = and excludes 
this function from the set ^k(r), the set becomes incomplete since it lacks a function of the 
translation symmetry of k . This problem was emphasized by ThoulessJil 



Since functions \I/k(r) are normalized, the denominator of \l/ko( r ) a l so turns into zero, 
z/(k ) = 0. Indeed, for k = k the cosine in Eq. (|7|) equals cos(k R mn ) = (— ) m+?1 ) and 
therefore z/(ko) = due to Eq. Function z/(k) is shown in Fig. [1] for a square lattice, 
a = b = \phx. One can see that z/(k) is positive, z/(k) > 0, inside the Brillouin zone and 
reaches its minima, z^(ko) = 0, in the corners of it. Both the numerator and denominator of 
\&k( r ) have the order of magnitude of |£| when £ = k — k — > 0. The leading term in the 
numerator of ^(r) depends on ip^ , the azimuth of £ , as it is shown below; cf. Eq. (|16|). 
On the contrary, the function y/ u(k), which is the denominator of ^k(r), is isotropic near 
ko for a square lattice. Therefore, the function ^/k(r) retains the dependence on (p^ even in 
the limit C, — > 0. Hence, it is singular at the point k = k , and the limit \I/k^ko does not 
exist. 

To find the form of \E , k(r) in the limit k — > ko , one can expand the numerator of \l/ko+c( r ) 
in £ and take into account Eq. (|4]). Then 

oo 

* k « + < (r) = \aJk)Y/i D { - r+n C ' R ^ Cm " (r) • (12) 

L " >• 'J mn=—oo 

If one introduces complex variables z = x + iy and Z mn = X m + iY n , the condition M) takes 



the form: 



oo 

£ (-) mn+m+n exp{--RL + ^^ m n} = 0. (13) 



^ — rrm ■ ^ 

mn=—oo 



Here and below complex conjugate variables are designated by bars. Since this equation 
is valid for arbitrary values of z and the sum converges exponentially, one can take the 
derivative over z: 

oo 

(-) m+n Z mn c mn (r) = . (14) 

mn=—oo 

If one takes advantage of the relation £ ■ R mn = (Z mn ( + Z mn ()/2, where ( = ( x + i( y , and 
plugs Eq. (0) into Eq. (|12|), the latter takes the form: 

. , oo 

^o+<(r) « [a ^( k)] i/ 2 C E (-) m+n Z mn c mn (r). (15) 



mn=—oo 



Since z/(k) 1//2 oc |£| for |£| — > 0, Eq. ( [T3|) can be rewritten as 



*k 0+ c(r) ~ e 



m *k (r) 



(16) 



where $k ( r ) is a regular function of r possessing the symmetry of the point k . 

Therefore, nonanalitic function ^ko+cl 1 ") factors near the singular point into the product 
of two functions, e~ llfi( and $k ( r ) • The first factor absorbs the nonanalytic dependence 
of \I/ko+c( r ) on C' whereas the second factor does not depend on £ and possesses the trans- 
lational symmetry of the point ko. Function $k (r), taken with an arbitrary phase factor, 
can be used as a Bloch function with the ko symmetry. Therefore, Bloch functions are 
defined for all k values. This statement concludes the proof of the completeness of the set 
^kM- The set W mn {r) obtained from it by the orthogonal transformation of Eq. (||) is also 
complete. 

There exists another way to construct the function <£>k (r). It is based on the properties 
of Bloch functions constructed from different sets of localized orbitals. Instead of coo(r), one 
can use the function c^(r) = zc o(r)/y/2 to generate the sets of orbitals c£i(r) and Bloch 
functions \l/^(r). The sets c mn (r) and Cmn(r) belong to the LLL. Both sets are overcomplete 
and can be expanded one in another. Consequently, functions \?k(r) and ^/^(r) can differ 
only in r independent phase factors for arbitrary value of k. Boon et al. have shown@ that 
k is a regular point for the set ^/^(r). Therefore, one can use ^/^(r) as a function $k (r). 
The shape of the function W mn (r) depends on the choice of the phase of $k ( r )- 

The properties of Bloch functions discussed in the previous paragraph can be also un- 
derstood from a more general point of view. It follows from Eqs. (|3|) and (|13D that functions 
Cmn (r) can be obtained as 



Therefore, functions c™(r) and c m „(r) belong to the same Landau level because of the 
existence of continuous group of magnetic translations. The same is true for the functions 
Cml(r) originating from the function c ^(r) oc z q Coo{r), where q is an arbitrary integer. Bloch 
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functions \E^(r), constructed from them, are completely determined by the translational 
symmetry up to the phase factors depending only on q and k. Therefore, all functions ^ (r) 
with the same value of the momentum k and different values of q coincide up to these phase 
factors. Singularities of functions ^^(r) ensure the existence of the Hall current J3 hence, 
they are present for each set of functions \l/^(r) and can not be eliminated. However, the 
number of singular points and their positions in the Brillouin zone change depending on g.@ 
It is known that the overcompleteness equation (^) is related to the properties of Jacobi 
§(u\t) functions.0'lll Eq. (|7|) for z/(k) takes a simple form when rewritten in terms of these 
functions. If one performs the summation over n in Eq. (0) using the Jacobi imaginary 
transformation of theta functions,i! i/(k) acquires the form: 



where k± = k x ± ik y . The zeros of ^s(u\r) can be found from the condition u = (M + |) + 
(N + |)t, where M and N are integers, and one immediately recovers that the corner of the 
Brillouin zone, ko = (±7r/a, ± 7r/6) , is a zero of ^(k). 

III. WANNIER FUNCTIONS IN THE INFINITE PLANE 

Equations (H) and (|ToD can be used to find Wannier functions w mn (r) localized in the 
infinite r plane. In the limit a, (3 — > oo, the sum in Eq. (|T0|) transforms into the integral 
over the Brillouin zone: 

„ ( v , f dk exp(ikR mn ) 
J(b.z.) ( 27r ) V^(k) 

The kernel Koo (m,n) is obviously independent of <p. 

The asymptotic behavior of K^m^n) for \m\, \n\ ^> 1 can be found analytically. It is 

determined by the behavior of the integrand near its pole, i.e., by the expansion of z/(k) 

near k . For a square lattice, a = b = a/27t, this expansion has the form 

u(ko + C) « 7« 2 C 2 /2 , (19) 



where 7 m 0.5814 is given by the series: 



7 = - n mn+m+n m 2 exp[~(m 2 + n 2 )}. (20) 



mn=—oo 



Substituting fll9|) into (|18|) results in the leading term of the expansion of K^^m, n) in -R m * : 

« h« 2 ,/| r « ^5=) = _y=L . (2D 



The next term of this expansion falls off as R^ n - An equation equivalent to K^m^n) oc 
(—) m+n /R mn was derived by Sen and Chitrall using a different procedure. 

The asymptotic behavior of w mn (r) for large |r — R mn | values follows from Eqs. @ and 

The right hand side of the equation 

00 

woo(r) = K^m.n) c mn (r) (22) 



mn=-oo 



for the function ^oo(r) includes the product of the kernel K^im, n), whose denominator 
depends on m and n slowly, and the factor c mn (r), which depends on m and n exponentially 
for a fixed value of r. One can use the asymptotic form of the kernel K^m, n), Eq. (pT|). 
and neglect the higher order corrections to it. If one substitutes the expansion 



iCn-^l-r-tR™,-!-)/ 



r 2 } 



into Eq. fl22|) , the first vanishes because of the Perelomov identity, Eq. @, and the second 
term results in 



7T7 r 3 

mn=-oo 



£ (-) m+n (R™ - r) c m „(r) . (23) 



The first factor in Eq. (|23"D falls off as r 2 , whereas the second is a periodic function of r. If 

one plugs r = R m v in the second factor and applies Eq. (|), this factor takes the form: 

00 

( _ )m '+„'+mV J2 (_)™+«+™» R mn Coo ( r ) . 
mn=-oo 

The sum is equal to zero because of symmetry arguments, therefore, all lattice sites, r = R mn , 
are zeros of the asymptotic expansion (|23|) of the function w mn (r). 
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It was shown by KohncZI that the rate of the falloff of an exponentially decaying Wannier 
function is determined by the distance of the singular point in the complex momentum plane 
from the real axis. In the problem of magnetic Wannier functions unavoidable singularities 
exist in the real (k x , k y ) plane. These singularities result in the power-law falloff of the 
Wannier functions. 

The asymptotic expansion of K^m, n), Eq. (plf), is accurate up to the values \m\, \n\ m 1, 
e.g., the deviation of K^l, 0) ~ —0.288 from its approximate value following from Eq. (|2~lf) 
is only about 2%. The coefficient i£oo(0, 0) ~ 1.241 is much larger than all the coefficients 
Kooim, n) with m,n ^ 0. One can subtract this large term and rewrite Eq. ( p2]) in the form 

oo 

w 00 (r) - c 00 (r) = (-) m+n A(m 2 ,n 2 )c mn (r), (24) 

mn=—oo 

The kernel A (m 2 , n 2 ) is numerically small and can be considered as a smooth function of m 
and n. This allows the application of the arguments employed when deriving Eq. fl2"3|). Then 
the leading term in the right hand side of Eq. ( f24"|) vanishes because of the overcompleteness 
condition of Eq. (Q). Therefore Woo(r) ~ Coo( r )> which means that the difference between the 
functions u>oo(r) an d c o(r) is expected to be small in the region where these functions are 
large. In the large r region, where these functions are small, iUoo( r ) dominates and Eq. ( |2"5D 
should be used. 

Numerical results support all of the above conclusions. The infinite-plane function i%](r) 
is shown in Fig. |2|. Fig. 0a provides a detailed comparison of the shapes of the function 
| w 00 (r) | 2 and the Gaussian function c 00 (r) 2 in the region of \x\, \y\ < 1.5 x \^2n. These 
functions are plotted along the x axis and along the diagonal, x = y, on the right and left 
hand sides of the figure, respectively. The function ^oo( r ) is re& l on these lines. One can see 
that u>oo( r ) 2 and coo(r) 2 are very close to each other in both directions. The function u7 o( r ) 2 
shows small anisotropy: it is slightly elongated in the x = y direction and squeezed in the 
x direction compared with Cg (r). It is a remarkable property of the function |w7 o( r )| 2 that 
it is very small in the points r = y/2n(l, 0) and v / 27r(l, 1). These points are the lattice 
sites adjacent to the origin. The suppression of |u> o(r)| 2 at the lattice sites of the two first 
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coordination spheres of the von Neumann lattice results in the squeezing of the central peak 
and a strong localization of |w o( r )| 2 - Consequently, the central region contributes about 
95% to the normalization integral. The behavior of woo{r) along the x axis is shown in 
Fig. ^]b over a wide region of x values. It is seen that woo(x, 0) oscillates and decreases with 
x. For x — > oo, the oscillation amplitudes decrease as x~ 2 and zeros of w 00 (x,0) approach 
multiples of the lattice period \/2ti in agreement with Eq. fl23l) . The exact u>oo(r) curve is 
almost indistinguishable from its asymptotic shape in the entire oscillatory region. 

Therefore, functions w mn (r) show a two-scale behavior. The function u>oo( r ) is large and 
shows only minor deviations from the Gaussian shape inside the central cell, but it is small, 
oscillates, and falls off according the r -2 law in the asymptotic region. 

The set w mn (r) is orthogonal and complete by arguments of Sec. ||. Therefore, any 
function which is invariant under the unitary transformations of the basis can be calculated 
in the Wannier representation. For instance, a straightforward calculation based on the 
representation of the functions w mn {r) through the kernel K^m^n) and coherent functions 
Cmn(r) shows that the function 

oo 

C 00 (r,r / )= ^2 w mn {r)w mn {r') (25) 

mn=— oo 

is equal to the expression 

CUr-r) = (l/27r)exp{-(r-r') 2 /4}exp{^ z- (r x r')}. (26) 



which can be also derived in the Landau representation. The last factor in Eq. (Pq) is 
gauge dependent. Therefore, the continuous symmetry of the magnetic translation group is 
recovered due to the completeness of the system w mn (r). 

IV. SOME RESULTS FOR TOROIDAL GEOMETRY 

Only few results can be obtained for Wannier functions W mn (r), Eq. (|8|), analytically. 
However, numerical methods allow one to study the dependence of W mn (r) on the plaquette 
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size (a, (3) and the twist 0, and to check the self-consistency of the procedure. A number of 
results was obtained by Ferrari.HH The singular behavior of ^k( r ) near k = k was studied 
in Sec. [TT]. 

In Fig. |] the square of the modulus of the function W / oo( r ) as well as the real and 
imaginary parts of Woo( r ) are shown for a square lattice with a = (3 = 3 and = 0. It is 
seen that |W / oo( r )| 2 is nearly isotropic and well localized in the area of about v^vr near the 
origin. The shape of this function is close to the data! for a triangular lattice. Re{W / oo( r )} is 
also rather isotropic and well localized, whereas Im{W / oo(r)} is small and highly anisotropic 
in the same region. The shape of the function Im{Woo( r )} can be understood if one takes 
into account that, because of Eqs. @ and @, the expansion of Im{Woo( r )} i n the powers 
of z starts with the term Im{Woo( r )} ^ Im{2; 4 }. All lower-power terms cancel for a square 
lattice. 

It is seen in Fig. |] that all three functions possess full symmetry of the square lattice. 
For an odd-odd plaquette this high-symmetry shape appears only for the twist = 0. When 
twist increases, the shape of |Woo( r )| 2 changes and becomes asymmetric. The changes are 
moderate near the maximum of the surface | WooC 1 ") | 2 but are considerably larger in its lower 
part. Dependence of the function W / oo( r ) on becomes discontinuous when one of the 
points determined by Eq. (g) passes through the corner of the Brillouin zone ko, since in 
the vicinity of k the function \E'k( r ) critically depends on the phase (p^ of £ = k — k as it 
follows from Eq. ([16]). 

More generally, the function |W m n(r)| 2 depends on disco ntiniously when i) both a and 
P are odd and <fi x = <p y = it, ii) both a and (3 are even and = 0, and iii) a is odd, (3 is 
even, and <f) x = tt and <p y = 0. 

Toroidal Wannier functions constructed by magnetic translations of infinite-plane func- 
tions w mn {r) have more stable shape than W mn {r). This approach will be discussed in more 
detail elsewhere. 

It was shown in Sec. [TJ that despite of the singular dependence of \l/k( r ) and W mn (r) on 
the twist in the vicinity of ko, these sets of functions remain complete at any values of 0. 
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The completeness can be checked by comparing the results obtained in the W mn (r) lattice 
basis with those in Landau function! basis. For example, the function 

C(r, r') = J2 W m n(r)W mn (r') , (27) 

mn 

which is similar to the function Coo(r, r ), Eq. (p5|) , does not depend on a specific choice of a 
complete basis set of the LLL. We have calculated the right hand side of Eq. fl27f ) for finite 
plaquettes using both W mn (r) and Landau functions as given in Ref. 2. Gauge independent 
parts of C(r, r') in both representations coincide with the accuracy of our computations. 
They were performed for the values of a and f3 up to 5. For the finite plaquettes the electron 
density p(r) = C(r, r) oscillates with r near its continuum limit = 1/2tt. The amplitude 
of oscillations decreases rapidly with increasing a and (3. The number of oscillations in x 
and y directions equals the number of fluxes af3 per a plaquette rather than the individual 
values of a and (3, which are the number of lattice sites in corresponding directions. This 
property was found analytically in the Landau representation by Sutherland.il It indicates 
disappearance of the pattern of the flux lattice and restoration of the symmetry of the 
underlying problem as it is expected since both periodic functions W mn (r) and Landau 
functions form complete sets. The completeness of the W mn {r) set is in agreement with the 
statements by Ferrari .1 



V. CONCLUSION 

The construction of infinite-plane localized magnetic Wannier-type functions was con- 
sidered a challenging problem for a long time. It was shown that these functions, if they 
do exist, are subject to rigid restrictions. Three properties of a set of localized magnetic 
orbitals are incompatible: completeness, orthogonality, and exponential falloff.0 Moreover, 
for systems supporting a Hall current magnetic orbitals should fall off no faster than r~ 2 .III 
The set of Gaussian coherent states on a von Neumann lattice violates these criteria because 
of the overcompleteness relation,0 which is a single linear constraint relating an infinite set 
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of orbitals. We believe that the complete set of orthogonal two-scale orbitals w mn (r) studied 
in this paper is the best compromise between the requirement of optimal localization and the 
inevitable restrictions on the degree of localization. Function u>oo( r ) centered at the origin 
i) is large and possesses nearly Gaussian shape inside the region of the size of about \phxl 
making a dominant contribution to the normalization integral, and ii) is small and falls off 
with a critical exponent r -2 outside this region. It is a striking property of the theory that 
the overcompleteness condition, Eq. emerges and plays a crucial role at all stages of the 
derivation and study of the Wannier orbitals w mn (r). 

We expect that w mn (r) orbitals form a convenient basis both for analytical and numerical 
calculations. 
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FIGURES 



FIG. 1. Normalization factor i/(k) for a square lattice. 

FIG. 2. (a) Dependence of the functions |u>oo(r)| 2 and coo(r) 2 on r in the directions of the 
principal axes x and y (right) and the diagonals x = ± y (left). Solid line - |u>oo( r )| 2 > dashed line 
- Gaussian function cqq{t) 2 . (b) Oscillatory dependence of the function wqq{x, 0) on x. The exact 
solution , asymptotic solution, and envelope function proportional to x~ 2 are shown in the region 
x > \f2nl by solid, dotted and dashed lines, respectively. 

FIG. 3. Shape of the function Woo(r) for a square lattice with a = (3 = 3 and = 0. (a) - 
|IFoo(r)| 2 , (b) - Re{VF 00 (r)}, (c) - Im{VF 00 (r)}. 
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